Different hierarchy of avalanches observed in Bak-Sneppen evolution model 
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The term of avalanche may originate from the phenom- 
ena occurred in nature. It is referred to as a sequential 
events which may cause devastating catastrophe. The 
phenomena of avalanches are ubiquitous in nature. The 
canonical example of avalanche is the mountain slide, 
during which great mass of snow and ice at a high al- 
titude slide down a mountain side, often carrying with it 
thousands of tons of rock, and sometimes destroy forests, 
houses, etc in its path [1]. Since avalanches occur ev- 
erywhere, from the ricepile, to the Himalayan sandpiles; 
from the river network, to the earthquake, starquakes and 
even solar flares; from the biology, to the economy, [2], 
etc, it is hence proposed [2] that avalanches may be the 
underlying mechanism of the formation of various geo- 
graphical structures and complex organisms, e.g., brains, 
etc. (It is now even proposed by Meng et al [3] that the 
formation of colorless gluon clusters may be attributed to 
avalanches intrigued by emission or absorption of gluons.) 
Prom this point of view, avalanches can be viewed as the 
immediate results of complex systems, and hence can be 
used as the theoretical justification for catastrophism. 
This is because if the real world is complex then the 
catastrophes are inevitable and unavoidable. 

Plenty of patterns provided by nature exhibit coher- 
ent macroscopic structures developed at various scales 
and do not exhibit elementary interconnections. They 
immediately suggest seeking a compact description of 
the spatio-temporal dynamics based on the relationship 
among macroscopic elements rather than lingering on 
their inner structure [4] . That is, one needs to condense 
information when dealing with complex systems. Maybe 
only this way is efficient and turns out successful. 

As known, avalanche is a kind of macroscopic phe- 
nomenon driven by local interactions. The size of 
avalanche, spatial and temporal as well, may be sensi- 
tive to the initial configuration, or more generally, the 
detailed dynamics of the system. However, the distri- 
bution of avalanches, i.e., Gutenberg- Ritcher law [5], or 
equivalently, power-law, does not depend on such kind of 
details due to the universality of complexity. Hence, in 
this sense avalanche study may be an appropriate tools 
in studying various complex phenomena. On the other 
hand, observation of a great variety of patterns, such 
as self-similar, fractal behavior in nature [6,7,8,9], 1// 
noise in quasar [10], river flow [11] and brain activity 
[12], and many natural and social phenomena, includ- 
ing earthquakes, economic activity and biological evo- 
lution suggests that these phenomena are signatures of 
spatiotemporal complexity and can be related via scal- 
ing relations to the fractal properties of the avalanches 
[13]. This suggestion means that the occurrence of these 
general, empirical phenomena may be attributed to the 
same underlying avalanche dynamics. Thus, one can see 
that study of avalanche is crucial in investigating the crit- 
ical features of complex systems. It can be even inferred 
that avalanche dynamics does provide much useful infor- 
mation for us to understand the general features of the 
ubiquitous complexity around us. 



Despite the fact that avalanche may provide insight 
into complexity, the definition of which can be vastly 
difii'erent for various systems, and the same sorts of sys- 
tems, even the same system. Let us recall some defini- 
tions of avalanche given before. In sandpile model [2], 
an avalanche is intrigued by adding a grain or several 
grains of sand into the system at some time and causing 
the topple of some sites, which may later on cause some 
other sites to topple. The avalanche is considered over 
when the heights of all the sites are less than the critical 
value, say, 4. In Bak-Sneppen model [14], several kinds of 
avalanche [13] are presented. For instance, /o-avalanche, 
G(s)-avalanche, forward avalanche, backward avalanche, 
etc. Though these kinds of definitions of avalanche may 
show various hierarchical structures they manifest the 
same underlying fractal feature of the system, i.e., self- 
organized criticality (SOC). Relating all these kinds of 
avalanches one can provide a general definition of the 
avalanche for Bak-Sneppen model: An avalanche corre- 
sponds to sequential mutations below certain threshold. 
One can see that this kind of definition can ensure the 
mutation events within a single avalanche are casually 
and spatially connected. In addition, with this defini- 
tion there exists a hierarchy of avalanches, each defined 
by their respective threshold. It is the hierarc:liic;al struc- 
ture of the avalanche, which exhibits the fractal geometry 
of the system and implies complexity. 

It can be inferred from the definition of avalanche that 
there always exists a triggering event which initiates an 
avalanche and whose effect, that is, causing the avalanche 
to spread later on within the system , will disappear 
at the end of the avalanche. Observation of avalanche 
through the triggering event, up to now, is based on the 
individual level, despite that the avalanche is a macro- 
scopic and global phenomenon of the system studied, in 
the laboratory, and in nature as well. In sandpile model, 
the triggering event is adding a grain or several grains of 
sand to some sites and causing them to topple. In Bak- 
Sneppen model, the corresponding triggering event of an 
avalanche is mutation of the extremal species causing the 
fitness [14] of the extremal site at the next time step 
less than a certain threshold. In the above two models 
triggering events are directly connected with the feature 
of individuals, e.g., the height of the site in the former, 
or, the fitness of the extremal site in the latter. It can 
be readily learned that the triggering events, whether 
in the laboratory or in nature, are not directly related 
to the global feature of the systems although avalanche 
can span the whole systems. Generally speaking, the be- 
havior of avalanches is observed through the features of 
individuals, instead of those of the whole system. How- 
ever, general features of the complex system may provide 
insight into knowing the tendency of the evolution of the 
system. Specifically, global features of a complex system 
may enable one to understand the critical behavior of 
the system. It implies that some characteristic quanti- 
ties, representing the corresponding global features, can 
be employed in describing complex systems. Further- 
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more, these quantities ought to be related to avalanche 
dynamics, and hence can be used to describe complexity 
emerged in a variety of complex systems. Apparently, 
our aim is to search for or define such kind of quantities 
and expect to observe new types of avalanches based on 
these quantities. Indeed, we obtain a new quantity which 
can be used to define a different hierarchy of avalanches 
in Bak-Sneppen model. We suggest that this quantity 
may be used as a criterion in determining the emergence 
of criticality. It will be shown later that this new type 
of avalanche still exhibits spatio-temporal complexity in 
another context. 

Bak-Sneppen model [14] is a very simple evolution 
model of biology. Despite the simplicity of the model 
itself, it can exhibit the skeleton of species evolution — 
punctuated equilibrium. In Bak-Sneppen model, each 
species is represented by a single fitness. The fitness may 
represent population of a whole species or living capa- 
bility of the species [15]. It is a vital quantity and the 
only one describing the model. No other additional quan- 
tities are employed in this oversimplified model. Thus, 
the fitness is the most important feature of species, and 
absolutely the most important of the model. So, when 
considering global feature of the ecosystem, one has to 
relate it to that of individuals, i.e., fitness. 

In this "toy" model (Bak-Sneppen model), random 
numbers, fi, chosen from a flat distribution, p{f), are 
assigned independently to each species located on a d- 
dimensional lattice of linear size L. At each time step, 
the extremal site, i.e., the species with the smallest ran- 
dom number, and its 2d nearest neighboring sites, arc 
assigned 2d + 1 new random numbers also chosen from 
p{f). This updating process continues indefinitely. After 
a long transient process the system reaches a statistically 
stationary state where the density of random numbers in 
the system vanishes for f < fc and is uniform above fc 
(the self-organized threshold). 

Introduce a new quantity for Bak-Sneppen model. De- 
fine the average fitness, denoted by /, as, 

i=l 

where fi is the fitness of the ith species. Here, we re- 
fer to / as the average fitness and a global quantity of 
the system. / may represent average population or av- 
erage living capability of the whole ecosystem. Large /, 
i.e., high average fitness, may imply the total population 
of the system is immense or its average living capabil- 
ity is great, and vice versa. Initial value of /, denoted 
by /(O) , can be easily obtained. As known, at the be- 
ginning of the evolution /j's are uniformly distributed 
between (0,1). Hence, for an infinite-size system, /(O) 
equals to 0.5. However, for a finite-size system /(O) fiuc- 
tuates slightly due to the finite size effect. It should be 
pointed out that /(O) does not reflect the correlation be- 
tween species. As the evolution goes on such correlation 
tends to be more distinctive. Denote f{s) the average 



fitness of the system at time step s in the evolution. In 
the s limit, i.e., s 3> L"^, f{s) may partly reflect informa- 
tion about correlation. As a global quantity, /(s) should 
include information concerning the interaction between 
species. Hence, it is natural to expect that / may be a 
good quantity in describing the feature of the system as 
a whole. 

Before introducing the different hierarchy of avalanches 
it is necessary and worthwhile to investigate features of 
the new quantity, /(s). Firstly, let us present some theo- 
retical analysis. Recall the definition of / one can see that 
A/(s) = f{s + 1) — /(s) approaches zero in the L ^ oo 
limit. An observer can hardly perceive the change in /(s) 
during such a short time period since it is vanishingly 
small. However, changes at each time step are accumu- 
lated to form a relatively distinctive change after a long 
time, which is perceivable for the observer. This long 
time period is required to be much larger than the sys- 
tem size, i.e., s L"^. In other words, f{s + Sq) — fi^o) 
may only be "noticed" when s ^ L'^ (sq denotes any 
initial time step). The variation of /(s) is small between 
two successive time steps, which differs from that of fit- 
ness of extremal site. The latter can be very large, say, 1. 
It should also be expected that there exists an increasing 
tendency of /(s) versus time s. This is because at each 
time step the least fitness is eliminated from the system 
so the general fitness of the whole system will tend to 
increase. And duc to the slow fluctuation of /(s) the 
increasing in /(s) behaves like a stepwise, i.e.. Devil's 
stepwise [2]. Hence, one may expect to observe such be- 
havior, i.e., pimctuatcd equilibrium [14], of /(s) in the 
evolution of Bak-Sneppen model. 

In order to show the feature of J{s) versus time s we 
perform simulations of Bak-Sneppen model. At each time 
step, in addition to the updating of the extremal sites, 
we also track the signals /(s). Fig. 1(a) presents the 
evolution of /(s) versus time s during some time period. 
This plot shows that /(s) varies slightly between two 
successive time steps but tends to increase in the long 
evolution process. 

Introduce another quantity, F{s), the gap of the av- 
erage fitness. The definition of F{s) is given as follows: 
Initial value of F{s) is equal to /(O). After s updates, a 
large F{s) > F{0) opens up. The current gap F(s) is the 
maximum of all F(s'), for all < .s' < s. Fig. 1(b) shows 
F(s) as a step-wising increasing function of s during the 
transient. Actually, the gap is an envelope function that 
tracks the increasing peaks in f{s). Indeed, punctuated 
equilibrium behavior appears in terms of f{s). 

By definition [14], the separate instances when the gap 
F{s) jumps to its next higher value are separated by 
avalanches. Avalanches correspond to plateaus in F{s) 
during which /(,s) < F{.s). A new avalanche is initiated 
each time the gap jumps and ends up when the gap jumps 
again. As the gap increases, the probability for the aver- 
age fitness, /(s), to fall below the gap increases also, and 
larger and larger avalanches typically occur. 

We can derive an exact gap equation of F{s), simi- 
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lar to that found in Ref. [16]. Suppose in the system 
the current gap is F{s). If F(s) is to be increased by 
AF, i.e., from F{s) to F{s) + AF, the average num- 
ber of avalanches needed is A^av = AFL'^/{fc — F{s)), 
where fc is the critical value of f{s). We can guaran- 
tee TVav > 1 by selecting AF > L~'^. In the large L 
limit, A^av can be arbitrarily large. Hence, in this limit, 
the average number of time steps required to increase 
the gap from F{s) to F{s) + AF is given by the inter- 
val As = (5)F(s)iVav = (S)F(,)AFL'^/(/e - F{s)), where 
(S')f(s) is the average size of avalanche of the plateaus 
in the gap function. /^From the law of large numbers 
the fluctuation of this interval around its average value 
vanishes. In the AF limit, As 0. Taking the 
continuum limit we can obtain the differential equation 
for F(s), 



dF(s) A - F(s) 



ds 



(2) 



Note this equation is exact. 

All SOC models, e.g., the BTW sandpile model [17], 
the earthquake models [18], or Bak-Sneppen model [14], 
exhibit self-organized criticality in terms of a power-law 
distribution of avalanche. It is natural to expect that we 
can observe SOC in terms of the hierarchical structure 
of /(s), which itself manifests complexity. It is simply 
another way to observe the same phenomenon by using 
such new quantity to define the avalanche, which can be 
observed in different ways. As known, the emergence of 
complexity is independent of the tools used to observe 
them provided that these tools are efficient and strong 
enough. Similar to those used in Refs. [13,19], we present 
the definition of /o-avalanche, where /o (0.5 < /o < 1.0) 
is only a parameter used to define the avalanche. Sup- 
pose at time step si, /(si) is larger than /p. If , at 
time step Si+l, /(si + 1) is less than /o, this initiates a 
creation-annihilation branching process. The avalanche 
still continues at time step s', if all the /(s) arc less 
than /o for 1 < s < s' — 1. And the avalanche stops, 
say, at time step si + S, when f{s\ + S) > /(si). In 
terms of this definition, the size of the avalanche is the 
number of time steps between subsequent punctuation 
of the barrier /o by the signal f{s). In the above ex- 
ample, the size of the avalanche is S'. It can be clearly 
seen from Fig. 1(a) that this definition guarantees the 
hierarchical structure of avalanches- larger avalanches 
consists of smaller avalanches. As Jo is lowered, big- 
ger avalanches are subdivided into smaller ones. Hence, 
the statistics of /o-avalanchc will inevitably have a cut- 
off if /o is not chosen to be fc- We can also define fc- 
avalanche. Nevertheless, /o-avalanche in the stationary 
state has the same scaling behavior as /c-avalanche pro- 
vided /o close to fc- We measure /o-avalanche distribu- 
tion for one-dimensional (ID) and two-dimensional (2D) 
Bak-Sneppen models. The simulation results are given 
in Fig. 2. The exponent r, defined by P{S) ~ S'"'^, 
is 1.800 for ID model and 1.725 for 2D model. An- 



other exponent, D, avalanche dimension [13], defined by 
n-cov S^^'^, where ricov is the number of sites covered 
by an avalanche, and d is the space dimension, is mea- 
sured. We find £)=2.45 for ID model and 3.10 for 2D 
model. 

Up to now, a question is still unsolved. It is the crit- 
ical value of /, fc- This may be troublesome if the sys- 
tem size is finite, but when we consider large L limit, it 
can be easily accomplished. Recall the evolution of Bak- 
Sneppen model, or the detailed research of the model of 
Ref. [13], the densities of sites with random numbers is 
uniform above G and vanishes below G when L ^ oo, 
where G is the gap of extremal site. One can readily 
obtain 



lim /(s) = lim 



1 + G(s) 



Prom Eq. (3) one can immediately obtain 

l + /c 



lim fc 



lim 



(3) 



(4) 



Hence, fc can be easily determined from Eq. (4). Using 
the results of fc provided by Refs. [13,20], one can ob- 
tain fc, 0.83351 for ID model and 0.66443 for 2D model. 
However, Eqs. (3) and (4) are not valid for a finite-size 
system, since one can not ensure the distribution of ran- 
dom numbers during a finite-size system is really uniform. 
Due to the fluctuation of /(s) it is extremely difficult to 
determine the exact critical value of / for a finite-size 
system. One may estimate fc for a finite-size system 
using the simulation. We find that this value weakly de- 
pends on the system size when it is enough large and 
fc will approach the corresponding value for infinite sys- 
tems. Fig. 3 shows the fluctuation of / for a ID model 
of size L = 200 near the critical state. We note, in this 
curve, / fluctuates slightly around some average value 
and does not tend to increase any more during a long 
time period. We may say that the system approaches its 
stationary state. In this sense, we suggest that / may be 
a good quantity in determining the emergence of critical- 
ity. That is, the great fluctuation of /min will not affect 
us to determine when we approach the critical state. We 
need only to know the feature of /. It is more reason- 
able and easily accepted since / is a global quantity and 
condenses information of the system and its components. 

Why we call the /o-avalanche a different hierarchy of 
avalanches? Firstly, this new type of avalanche is defined 
on the global level, in terms of the new global quantity, /. 
The background of this definition is different from those 
used before. This type of avalanche reflects the fractal 
geometry in terms of the global feature. Secondly, one 
can notice that the exponents r of avalanche distribu- 
tion obtained in our simulations are different from those 
found in Ref. [13]. Prom this point of view, one can 
conclude that this type of avalanche is different from any 
one observed before. 
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SOC was suggested by Bak et al. to be the "finger- 
prints" of a large variety of complex system, which is rep- 
resented by a scale- free line on a double logarithm plot. 
In order to know the criticality of a system one needs 
to know when the system reaches the stable stationary 
state where the phase transition occurs. It is extremely 
difficult and nearly impossible for one to know when 
a system in nature approaches its critical state. One 
has to study the ubiquitous fractal geometrical structure 
carved by avalanches through thousands of millions of 
years. However, in laboratory experiments and computer 
simulations, one needs a criterion to determine whether 
stationary state approaches, even, reaches, since statis- 
tics of avalanches may only be done under critical state 
of the system. Given Bak-Sneppen model, when the ex- 
tremal signal, /mill, approaches the self-organized thresh- 
old, /c, the ecosystem reaches its stationary state. How- 
ever, /min itself fluctuates greatly time to time, which 
brings a big problem in determining the appearance of 
criticality. Thus, we provide a new quantity, /, for a can- 
didate in judging the emergence of criticality. As shown, 
/ is relatively stable in a short time period. Hence, when 
/ does not tend to increase any more, one may say that 
the system approaches its stationary state. And, one can 
observe criticality in a quite long time period. 

In conclusion, a different hierarchy of avalanches is ob- 
served in Bak-Sneppen model. A new quantity, /, is 
presented and suggested by us to be a possible candidate 
in determining the emergence of criticality. An exact gap 
equation and simulation results are also given. 
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Figure Captions 

Fig. 1: (a) The variation of / versus time during a time 
period for a one-dimensional Bak-Sneppen model of size 
L = 200. This shows the hierarchical structure of /. (b) 
Punctuated equilibrium of / for a one-dimensional Bak- 
Sneppen model of size L = 200. We track the increasing 
signal of /(s), i.e., F(s). 

Fig. 2: Distributions of /o-avalanche for (a) one- 
dimensional Bak-Sneppen model with size L = 200, 
/o = 0.821 and slope=-1.800; and (b) two-dimensional 
Bak-Sneppen model with size L — 20, /o — 0.648 and 
slope=-1.725. 

Fig. 3: The fluctuation of / around the critical state 
for a one-dimensional Bak-Sneppen model of size L = 
200. 
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